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Within a recently proposed variational approach it has been shown that, 

in Lai-xAxMnO-i perovskites with < x < 0.5, near the metal-insulator 
transition, the combined effect of the magnetic and electron-phonon inter- 
actions pushes the system toward a regime of two coexisting phases: a low 
electron density one made by itinerant large polarons forming ferromagnetic 
domains and a high electron density one made by localized small polarons 
giving rise to paramagnetic or antiferromagnetic domains depending on tem- 
perature. Employing the above-mentioned variational scheme, in this paper 
spectral and optical properties of manganites are derived for x = 0.3 at dif- 
ferent temperatures. It is found that the phase separation regime induces a 
robust pseudogap in the excitation spectrum of the system. Then the conduc- 
tivity spectra are characterized by a transfer of spectral weight from high to 
low energies, as the temperature T decreases. In the metallic ferromagnetic 
phase, at low T two types of infrared absorption come out: a Drude term 
and a broad absorption band due respectively to the coherent and incoherent 
motion of large polarons. The obtained results turn out in good agreement 
with experiments. 
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In the last years the perovskite oxides Lai^^AxMnO^ {A represents a divalent alkali ele- 
ment such as Sr or Ca) have become one of the main research areas of the condensed matter 
community because of the colossal magnetoresistance effect exhibited for 0.2 < a; < O.S.Bi 
The study of these materials started in the 1950's showing the strong correlation between 
magnetization and resistivity.i Since in the electronically active Mn 3d orbitals the mean 
number of d electrons per Mn is 4 — x, three electrons go into the t2g core states and the 
remaining 1 — x electrons occupy the outer shell eg orbitals. The ferromagnetic phase was 
explained by introducing the double exchange mechanismSi that takes into account the com- 
bined effect of the Cg electron hopping between nearest neighbor sites and the very strong 
Hund's exchange with the localized t2g electron spins. In the 1990's the discovery of the 
colossal magnetoresistance phenomenon has aroused a renewed interest for these compounds. 
In order to explain the experimental data, it has been suggested that, in addition to the 
double-exchange term, a strong interaction between electrons and lattice distortions plays 
a non negligible role.i Next the relevance of the Jahn- Teller polaron formation has been 
confirmed experimentally by the giant isotope shift of the Curie temperature,^ by measure- 
ments of the lattice distortions by EXAFsi and by frequency shifts of the internal phonon 
modes.lHl Furthermore both pseudogap features and conductivity spectra of these compounds 
have been discussed in terms of a strong coupling to lattice distortions.li2FEl Nowadays the 
role of the polaron is widely recognized.0 

Recent studies have showed the presence of strong tendencies toward the phase separation 
in manganites.i0 At low temperatures for compositions near the metal-insulator transition 
the ferromagnetic and antiferromagnetic phases coexist while, close to the paramagnetic to 
ferromagnetic transition temperature, the phase separation occurs between paramagnetic 
and ferromagnetic domains. From a theoretical point of view, the coexistence between 
hole-poor and hole-rich phases has been discussed by using exact numerical approaches on 
small lattices assuming classical Jahn- Teller phonons.Blii The pseudogap formation^ and 
the optical propertieJ^3ill also have been studied neglecting the quantum nature of phonons. 
One of the drawbacks of this approximation in the optical properties is, for example, that 



at low temperatures in the metallic ferromagnetic phase a narrow Drude peak cannot be 
obtained. 

In a recent paper@, some of us have shown, on the basis of a variational approach, that 
the quantum character of the Jahn- Teller phonons and the polaron formation can be impor- 
tant to explain the experimentally observed tendency of manganites to form inhomogeneous 
magnetic structures near the phase boundaries. We employ the scheme proposed in that 
paper to deduce spectral and optical properties of manganites for x = 0.3 at different tem- 
peratures. It is interesting to study the properties of these compounds at this doping since 
in the range of concentrations around x = 0.3 the colossal magnetoresistance effect is very 
pronounced in many manganese oxides.0 In this paper we find that the phase separation 
regime induces a robust pseudogap in the excitation spectrum of the system. The pseu- 
dogap features show a direct relationship with the metal-insulator transition and with the 
crossover from the coherent large polaron dynamics to the incoherent small polaron one. 
Our results turn out compatible with experimental findings.0'0 Next, by using the formal- 
ism of generalized Matsubara Green's functionsiHH^, we determine the scattering rate of 
charges affected by the interaction with the lattice and the spin fluctuations. The single 
phonon emission and absorption represent the main mechanism of scattering. However, the 
damping due to spin fluctuations is effective in the energy range around the chemical po- 
tential fi. The scattering rate turns out fundamental to derive the very interesting optical 
properties of the system. With decreasing T, our spectra are characterized by a transfer of 
spectral weight from high to low energies filling up the low frequency optical gap present in 
the high-temperature phase. Indeed at high T the infrared absorption is due to the incoher- 
ent small polaron dynamics. Instead, at low temperatures, in the ferromagnetic phase the 
system shows two types of optical response: a Drude term and a broad absorption band due 
respectively to the coherent and incoherent motion of large polarons. The results obtained 
are consistent with experimental data.BElii 

In section 1 the model and the variational approach are reviewed; in section 2 the spec- 
tral properties are deduced; in section 3 the damping of the particle motion is calculated; 



in section 4 the optical properties are discussed. In appendix details of the calculations 
employed to derive the optical properties are reported. 



I. THE MODEL AND THE VARIATIONAL APPROACH 



The model takes into account the double-exchange mechanism, the coupling of the eg 
electrons to lattice distortions and the super-exchange interaction between neighboring local- 
ized electrons. The coupling to longitudinal optical phonons arises from the Jahn- Teller 
effect that splits the Cg double degeneracy.!! Adopting the single orbital approximation 
(reasonable for x < 0.5), the Hamiltonian reads 



Here t is the transfer integral between nearest neighbor (nn) sites < i,j > for electrons 



on nn sites and the conduction electron, S is the spin of the t2g core states (S* = 3/2), c] (q) 
creates (destroys) an electron with spin parallel to the ionic spin at the i-th site. The first 
term of the Hamiltonian describes the double-exchange mechanism in the limit where the 
intra-atomic exchange integral J is far larger than the transfer integral t. Furthermore in 
eq.(|l]) ujq denotes the frequency of the local optical phonon mode, a| (oj) is the creation 
(annihilation) phonon operator at the site i, the dimensionless parameter g indicates the 
strength of the electron-phonon interaction in the Holstein modefi, e represents the an- 
tiferromagnetic super-exchange coupling between two nn t2g spins and /i is the chemical 
potential. 

The hopping of electrons is supposed to take place between the equivalent nn sites of 
a simple cubic lattice separated by the distance \n — n'\ = a. The units are such that the 
Planck constant h = 1, the Boltzmann constant kB=l and the lattice parameter a=l. 




(1) 



occupying Cg orbitals, Sq-' is the total spin of the subsystem consisting of two localized spins 
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Following the recently proposed variational schemeii, we perform two successive canon- 
ical transformations to treat the electron-phonon interaction variationally. 
The first is the variational Lang-Firsov unitary transformationiilii 



Ui = exp 



(2) 



where / and A are variational parameters. The quantity / represents the strength of the 
coupling between an electron and the phonon displacement on the same site, hence it mea- 
sures the degree of the polaronic effect. On the other hand, the site-independent parameter 
A denotes a displacement field describing static lattice distortions that are not influenced 
by the instantaneous position of the electrons. 

The second canonical transformation is Bogoliubov-typell 



U2 = exp 



—a ^ (alal — a^a, 



(3) 



where a is a variational parameter. It introduces correlations between the emission of 
successive virtual phonons by the electrons and it is responsible of the phonon frequency 
renormalization. 

The transformed Hamiltonian H = U2UiHUi^U2^ becomes 

H = ^i,A^j + '^0 1] ci\a"i + Nuq sinh^ (2a) + Nuog^A^ 



<i,3> 



+ujq sinh (2a) cosh (2a) ^ {ci\a\ + ajOj j — oooge 

i 

<i,j> i 

where the operator reads 

= —t ( „ — /— I xjx. 



(4) 



2S + 1 



(5) 



and Xi is the phononic operator 



Xi = exp gfe (a^ - aj) 



(6) 



Furthermore in the Hamiltonian (^ we specify the renormalized phonon frequency ujq 
ciJo cosh(4a), the number of lattice sites N, the phonon operator Ci 



C, = cuo(7(l-/)e2-(a, + al) (7) 

and the quantity rj 

rj = uog'f (/ - 2) + 2uog' (/ - 1) A (8) 

that measures the electronic band shift due to the electron-phonon interaction. 

To obtain the free energy in a variational scheme, we introduce a test Hamiltonian 
characterized by electron, phonon and spin degrees of freedom non mutually interacting 

Htest = -teff +^oYl ^l^i + siuh^ (2a) + Nuog^A^ 

<i,j> i 

-gslJ'sYKff Si + {j] - n)Yc\ci. (9) 

i i 

The quantity te// denotes the effective transfer integral 

'--{(^)>-- (-) 

where the symbol <> indicates a thermal average and the quantity St is 

St = g^pe-'^ {2N, + 1) (11) 

with A^o the average number of phonons with frequency ujq. In the test Hamiltonian (^ gs 
is the dimensionless electron spin factor {gs — 2), fis is the Bohr magneton and he// is the 
effective molecular magnetic field in a cell containing two neighboring sites. The magnetic 
field is induced by a partial orientation of ionic spins around this cell and is determined 
by the variational approach.0 We notice that in eq.(|TD|) the factor e^"^^ controls the band 
renormalization due to the polaron formation. 

To obtain the variational free energy of the system, we employ the Bogoliubov inequality 

F < Ftest +{H- Htest)t (12) 

where Of indicates a thermodynamic average made using the test Hamiltonian. This 
approach allows to treat the local spin dynamics in a variational mean field theory. The 
free energy per site becomes 



^ = ft, + T log (l - e-^^^) + uj, sinh^ (2a) + u,g\l - ffp" - T log us 

±^ZS^ml + TXms (13) 

where the electron free energy reads 

ftest = (^) E a^ik + ^ (^) E [^k log nk + (1 - nk) log (1 - nk)] (14) 

with n\f^ = Up (^k) the Fermi distribution function. We get ^k = — where ^k = + ^ 
is the renormalized electronic band and ^k the band dispersion 

^k = — 2te//[cos(fc^) + cos{ky) + cos{kz)\. (15) 

Furthermore in eq.(0) j3 is the inverse of the temperature, vs is the partition function of 
the localized spins, the top and bottom sign of e hold, respectively, for the ferromagnetic 
and antiferromagnetic solutions of the localized spins, Z indicates the number of nearest 
neighbors, represents the magnetization of the localized spins and A is a dimensionless 
variational parameter proportional to the effective magnetic field. 

In order to simplify the calculations, we consider a semicircular electronic density of 
states 

g{i) = (^) e{W-\i-r^ + ^^\)^W^-{^-^ + ^f (16) 

where W = Zteff is the renormalized band half-width and 6{x) is the Heaviside function. 
Actually g{^) represents a simple approximate expression for the exact density of states and 
it is generally used for a 3-D lattice.!!!! Therefore the electron free energy becomes 

/oo 
-oo 

/oo 
d^) {npiO lognp(0 + [1 - npiO] log [1 - n^(0]} ■ (17) 
-oo 

As shown in ref.@, in the intermediate electron-phonon coupling regime, the free energy 
(P!B|) gives rise to a region of coexisting phases characterized by different electron densities 
pi and p2 but the same value of chemical potential fi. The quantities pi and p2 depend on 
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the temperature and correspond to homogeneous phases of the system constituted by large 
and small polarons, respectively. Hence, near the metal-insulator transition (see Fig. la), 
the system segregates in ferromagnetic and antiferromagnetic or paramagnetic domains of 
large and small polarons. We observe that the quantity x denotes the hole concentration, 
that is X = 1 — p, where p is the electron concentration. 

In the region of coexistence, the fractions of volume, Vi/V and V2/V, filled with density pi 
and p2, are determined by the two conditions: Vi/V + V2/V = 1 and (yi/V)pi + (y2/V)p2 = 
p. In Fig. lb we report the fractions of volume (Vi/V) j:^^^^ and (V2/V)p^j,^ for x=0.3 as a 
function of the temperature (the dotted line in Fig. la allows to recognize the different phases 
of the system at this doping). 

In the regime of phases coexisting,0'ii it is possible to calculate any property B of the 
system by means of the linear combination of the properties Bi and B2 characteristic of the 
single phases assuming the respective volume fractions (Vi/V) and (V2/V) as weights 



Hence we assume that in this regime the properties of the system are independent of mor- 
phology of coexisting domains. 



In this section we calculate the spectral properties of the system. 

Performing the two canonical transformations (|^,§) and exploiting the cyclic properties 
of the trace, the electron Matsubara Green's function becomesH 




(18) 



II. SPECTRAL PROPERTIES 



G (k, r) = - (1) e^H^^^^s) (T^g^ (^) X,(r)4x}) 



(19) 



where (r) = e'^^Cie~'^^ and Xi (r) = e'^^ Xie~'^^ , with H given by eq.(^). 

The correlation function in (|19|) can be disentangled into electronic and phononic terms 
by using the test Hamiltonian (^, hence 
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Q (k, r) = - (1) e^'^-(^-^^) (T.Q (r) c]),(r.X, (r) X]), 



(20) 



where now we have Cj (r) = e^^^CiC '^^^ and Xj (r) = e'^^'XjC 
The Green's function in Matsubara frequencies Un reads 

e"""" (^) E (;^) (ki,z^^) |^^e(^'^--^--)^{e^^°^M^°(^^^)] - 1} (21) 

where Q^^^ (k, iun) is the Green's function of non interacting particles and s is 



s = 2/Ve-^" [iVo(iVo + l)]^. 



(22) 



We notice that two physically distinct terms appear in the Eq.(plD, the coherent and the 
incoherent one.00 The first derives from the coherent motion of electrons and their sur- 
rounding phonon cloud (this process is called diagonal transition) . The renormalization fac- 
tor e""^"^ appearing in this term plays the role of the residue Z at the pole in a Fermi-liquid 
description.B The second contribution describes the possibility of changing the number of 
phonons in the phonon cloud during the electron hopping (this process is known as non di- 
agonal transition) . This is confirmed by the presence of the sum over all electron momenta 
which reflects the fact that the electron momentum is not conserved separately from the 
phonon one. 

In the second term of the expression (0), we can perform the expansion 



exp < s cosh 



UJo T 



2^/;(s) cosh 
1=1 



UJo T 



+ Io{s) 



(23) 



into multiphonon contributions, where the symbols //(s) denote the modified Bessel func- 
tions. 

Making the analytic continuation iujn w + iS, one obtains the retarted Green's function 
G'ret(k, cu) and the spectral function 



A{k, u) = -2'^Gret{^, to) = 2ne-^^5 {u - ^k) + 
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27re"^ 2^ Ii{s)e^^[l — npi^uj — lLDQ)]g{uJ — Iuq) + 

/(^O)-oo 

27re~^ 2^ /i(s)e^~nir(co' + /c<}o)fi'(t^ + too) + 
27re-^-[/o(s)-l](7H. (24) 

The first term represents the purely polaronic band contribution and shows a deha behavior. 
The remaining terms provide the incoherent contribution and spread over a wide energy 
range. Thus we get the renormahzed density of states 

iVM = ^(^)E^(k,a;), (25) 

where the spectral function 74(k, cj) is given by eg. (p^) . We note that the sum rule 
^A(k, cu) = 1 is satisfied and N{uj) is normahzed to unity. 
So we can evaluate the renormalized electron distribution function 

duj 



/ + 00 f-j^^j 
—A{k,uj)nF{uj) = e-^^n^(ek) + 

,1 /n\ J— oo 



/{7^0)-oo 



l{^0)-oo -^""^ 

/oo 
dmWiS). (26) 
-oo 

At zero temperature it is found n'"^"(/i— 5)— n™"(/x+(5) = e""^^, so the factor e~'^^ determines 
the jump in the Fermi distribution function.0 For low temperatures and x > 0.25, e~'^^ is 
about 0.9, therefore the spectral properties indicate a polaronic Fermi-liquid state in the 
metallic ferromagnetic phase. 

On the other hand, for the small polaron excitations and high temperatures, we expect 
that the spectral properties are dominated only by the incoherent term. It is well knownii 
that in this limit the spectral function (|2^ ) consists of a sum of delta functions on the points 
nujo- This is caused by the choice of a dispersionless phonon frequency. To overcome this 
drawback, there are two ways of proceeding: the first is to introduce a phononic dispersion 
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relation^, whereas the second is to make a high temperature expansion.Ei In this paper 
concentrate on the latter approach performing in eq.(^T]) the following approximatio ni 



we 



e exp < s cosh 



exp[—sT) exp I — r 



(27) 



where s is 



4sinh(^) 



(28) 



This expansion is valid as long as (^^"j (^^) < hence T > O.Stoo- In this limit the 
Green's function becomes 



g (k, tUn) = e-^-6;(°) (k, tUn) - e-^- i^-j (ki, ^^n) + 



where W{z) is 



ki 



2Vi - ?ki - s) 



+ 



2\s 



TT JO 



W^(z) = e- 

In the small polaron case, we get ~ —g'^u^o — A*. The incoherent term yields 



(29) 



(30) 



A{k,uj)c^A{iu) 



l-np{-g ujQ- IJ, 
V s 



ln(3 
— exp 
s 



As 



+ 



f I 2 I i -\2 

- — [uj + g ujQ + ii + s) 
As 



(31) 



We stress that this spectral function represents the envelope of the delta functions obtained 
by eg. (p^) for small polaron excitations in the limit of high temperatures. Furthermore, 
using the spectral function (|3T|), we can derive the renormalized density of states 



N{uo) 



2ti 



(32) 
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that turns out normalized to unity. 

In Fig. 2 we report the spectral function for two different temperatures at x = 0.3. At 
T = O.OStuo, the spectral function is deduced by eg. p^) and the case .^k = is considered. 
It is apparent that in the large polaron ferromagnetic phase the weight of the coherent 
term is prevalent. For T = 0.470^0, in the paramagnetic phase, the spectral function is 
provided by eq.(|3l|) and it is determined by the incoherent dynamics of the small polaron 
excitations. It is made up of two bands peaked approximatively around —q'^ujq and g'^uo, 
whose different heights are due to the fixed doping. At high temperatures the spectral 
function is characterized by a large depression in correspondence of fi. With decreasing 
temperature and fixed doping, a transfer of spectral weight occurs crossing the regime of 
coexisting phases. To stress this important point, in Fig. 3 we present the renormalized 
density of states for x = 0.3 at different temperatures. In the large polaron ferromagnetic 
phase, N{u!) is given by eg. (pS]) , in the small polaron paramagnetic phase by eg. (|3T|) , in the 
regime of coexisting phases by the two preceding densities by means of the eg. (p!8|) . We stress 
that Tc denotes the ferromagnetic transition temperature which coincides with the beginning 
of formation of large polaron domains. For x = 0.3, is about 0.460^0. For low temperatures 
incoherent processes of phonon creation and annihilation lead to a depression in the density 
of states around the chemical potential /i. When the temperature increases, the effect of 
band renormalization driven by e""^^ is stronger than this decrease of electronic states at 
fi. Since the density of states is normalized, the net effect is a rise of states at /i (compare 
graphs at T = 0.28Tc and T = 0.65Tc). However, when the phase separation regime takes 
place (starting from T = 0.65Tc), the formation of paramagnetic insulating domains in the 
metallic background causes a robust pseudogap in the excitation spectrum. The subtraction 
of levels below /i is strongly temperature dependent and, as we see in the inset of Fig.3, it 
finishes at Tc. Therefore the pseudogap features show a direct relationship with the metal- 
insulator transition and the coherent-incoherent crossover of charge dynamics. Above Tc 
the density of states at fi is low and with increasing temperature it is nearly constant. It is 
important to note that these results are compatible with photoemission experiments. Indeed 

12 



the density of states at /i in the inset of Fig. 3 bears a strong resemblance with the spectral 
weights reported in refdl§0. 

Finally it is interesting to study the pseudogap formation fixing the temperature and 
varying the hole concentration x (see Fig. 4). At T = 0.39^70, the system is in the insulating 
paramagnetic phase from x = to x = 0.17, in the regime of the insulating paramagnetic 
and metallic ferromagnetic coexisting phases from x = 0.17 to x = 0.33, in the metallic 
ferromagnetic phase from x = 0.33. Increasing x, holes are added to the system and fer- 
romagnetic large polaron domains are created yielding pseudogap features. As soon as the 
phase separation regime is set up at x = 0.17 (inset of Fig. 4), electronic levels are occupied 
in the energy range around the chemical potential /x. Thus, the charge carriers prefer to be 
located in the metallic ferromagnetic domains in order to increase the kinetic energy gain. 



III. DAMPING 

In this section we deal with the polaron self-energy. This quantity allows to determine 
the scattering rate of particles affected by lattice distortions and spin fluctuations. The 
latter, in turn, will play an essential role in the infrared absorption calculations. 

Retaining only the dominant autocorrelation terms at the second step of iteration&0, 
we can derive the self-energy 

S(2) (k, lun) = e^+ (1) [^^ ^° (ki, ^^m) [ dre^^^-^^iMr) + Mr)] (33) 
where £k is given by eg. ([T5|) , /i(t) is 

with 

z = 4g'fe-'"[No{No + l)]"^ (35) 

and /2(t) is 
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/2(r) = 2gW^e^'^{l - fY[No{No + 1)]^ cosh 



UJQ T 



(36) 



Making the analytic continuation iujn —>■ w + i6, we can calculate the scattering rate 



r(k) = r(k, = a) = -2^JsS (k, ^ = a) • 



(37) 



We realize that the dependence of r(k) on k is due only to ^k- Hence, performing the 
expansion (|23D into the series of multiphonon processes, we can express the scattering rate 
in the following way 



r(k)=r(a) = ri_phon(ek) + r 



multi—phon 



(^k) + ^Spin-FluctiCk) 



(38) 



where Ti_phon is the contribution due to single phonon processes only 



r,_,..„(a)=2ZtV^^-^(^|±l/^^ 



Ii{z) sinh ( ^ I fi'i,/=i(^k) + 



(39) 



^muiti-phon represents the scattering rate by multiphonon processes 



^ multi—phon{^]i) IZt 6 

and Tspin-Fiuct denotes the damping term by spin fluctuations 



(40) 



Spin—Fluct 



25 + 1 



^ ^0 + 1/2 ^ 
2S + 1. 



B{i^.). (41) 



Here gi^i reads 

^i.K^k) = [A^o(too) + ni.(a + /cuo)] 5(a + l^o) + [NoilOJo] + 1 - np^^i, - too)] 5(a - too) 

(42) 



and g2 



^72(ek) = [A^o + npiii, + cuo)] 5(ek + c^o) + [A^o + 1 - ^^(ek - <:^o)] S(ek - cuo). (43) 



The function i?(x) is defined as 
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B{x) = 27ig{x) (44) 

where g{x) is the density of states (p^ . 

We stress that the decomposition of the scattering rate in three distinct terms has been 
introduced in order to simphfy the analysis of our results. Furthermore we observe that 
the scattering rate for small polaron excitations can be correctly calculated only within a 
self-consistent treatment.00 This can be carried out substituting in eq.(p3D Q^'^^ for a new 
Green's function Q. Using the Lehmann representation we can write 

g{k,tuJn)= — . (45) 



oo 2tT iuJn — UJ 



where the spectral function A is assumed to be 

So, employing this self-consistent procedure, we get the scattering rate r'^*'^)(^k) 

= ri_p/io„(^k) + '^[nulti-phoni^^ + '^'^Spin-Flucti^^ i^"^) 

where T^i%honi ^^^Im-phon and T^spin-Fiuct are calculated by eqs (|3|),(0),(||), respectively, 
substituting the function B{x) for B^^''\x) 

ki 

Since the spectral function A (^61) depends on F, the eq.(^7D allows to evaluate self- 
consistently the damping due to polaron formation and spin fluctuations. Clearly when 
in the right-hand side of eq.(p7D the limit F — > is made, the equation (|38D is obtained. 



B'^''\x) = (1) $:i(ki,x) = J^^d^M^Mi^ux). (48) 



In Fig. 5a we report the scattering rate F of large polaron excitations for different tem- 
peratures. Because of the augmented number of phonons and the enhanced role of spin 
fluctuations with rising temperature, the quantity F increases. At low temperatures we no- 
tice that the scattering rate is zero for energies within ujo of the chemical potential fi. This 
means that at these temperatures the spin fluctuations are too feeble to scatter electrons 
and the main mechanism of energy loss is the emission or the absorption of phonons with 
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frequency uq (one phonon processes are prevalent). Thus the behavior of the scattering rate 
is explained by the fact that the quasi-particle excitations within ljq of the Fermi energy 
cannot lose energy because of the Fermi statistics.^ 

In Fig. 5b we concentrate on the components of the scattering rate at a fixed temperature 
(T = 0.65Tc). It is confirmed that also at higher temperatures the single phonon emission 
and absorption give an important contribution to the quantity F. Nevertheless we stress 
that the spin fluctuations are effective to enhance the scattering rate just in the important 
energy range around the chemical potential /i. Finally the multi-phonon contribution is 
dominated by two-phonon processes and it is negligible. 

We stress that for the small polaron excitations the quantity F decreases when T increases 
resembling the behavior of tg//- We note, however, that, for T > O.Sujq, F is always larger 
than teff, making clear that the electronic states lose their individual characteristics and 
the electron motion is predominantly a diffusive process. 

The introduction of the damping allows to improve the approximations of calculation for 
the spectral properties: indeed we can determine the spectral function 

A{k,uj) = e-^^A{k,uj) + 

2e-^- ^/Ks) sinh (^^^ (1) [^o(too) + np{uj + luo)] A{k,,u + luo) + 

e-^-[Jo(.)-l](^)Ei(ki,c.). (49) 

In the limit F ^ 0, the result of eq. (|2^) is recovered. We have checked that the quantity 
F does not change the spectral properties in a considerable manner, although it allows to 
eliminate the delta behavior in the expression (^4|). 

On the other hand, as mentioned, the damping is essential for determining the optical 
properties of the system. 
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IV. OPTICAL PROPERTIES 



In this section we focus our attention on the optical properties. 

To determine the hnear response to an external field of frequency uo, we derive the 
conductivity tensor cTq ,g by means of the Kubo formula 



<^o.A^) = (^) < t,/^ > + K%^) (50) 

where T^,^^ is the kinetic energy and n^'^^(u;) is the retarted current-current correlation 
function. The transport of electrons will be supposed to take place between the equivalent 
nn sites of the cubic lattice, hence the tensor aa,[3 is assumed to be diagonal with mutually 
equal elements aa,a- To calculate the infrared absorption, we need only the real part of the 
conductivity 

^aa,a{^) = ■ . (51) 

Therefore our problem reduces to evaluate the current-current correlation function that in 
Matsubara frequencies is defined as 

n«,a(^u;n) = (-^) ^^^-"-(T.£(r)j„(0)) (52) 

where the current operator suitable for the Hamiltonian is (P 

Ja = lie }^ \^ — ^^-^^ j cl^s&Ci (53) 

and Jair) = e^'^j^^e"^'^ . Performing the two canonical transformations (0,^ and making the 
decoupling of the correlation function (^) in the electron, phonon and spin terms through 
the introduction of Htest @, we get 

r/3 . 



where the function A 



^'(-^)EE(^-^') / e^^"^A(^^',5,5')(T4(T)QH-..(r)4^,,,c,)* (54) 

^ i,5 i',S' ''^ 



A{i, i', 6, 6') = $ (z, i', 6, 6') T{i, i', 6, 6') (55) 
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denotes the product of the phonon correlation function $ 



times the spin correlation function T 

/ /^M+5a + l/2\ / ^^'+'5'a/ ^ 



We calculate the electron correlation function 



2S+1 



(56) 



(57) 



{Trcl{T)ci+s&{'r)4+s'&Ci')t 



k,ki 



and we separate A into two terms 



hence 



^{i,i',S, 6')T{i,i',6, 6') 



g^,^+,a ^^^^^ 
2S+1 . 



+ 



$ (i, t', 5, 5') T(z, z', 5, 6') - \{XjXi+sa) 



2 / p^'^+'" + l/2 ^ 
U 25+1 , 



(58) 



(59) 



A{i, i' , S, S') 



'2St 



' So + 1/2 ' 
. 25+1 , 



6')T{i,i',6, 6')-e 



_25, / / go + 1/2 ^ 
\\ 25 + 1 



• (60) 



To derive the optical properties, the role of the damping F of the particle motion is 
fundamental. Furthermore the defect scattering rate Fq is added to the quantity F calculated 
in the previous section. This scattering rate Fq is assumed to be a small perturbation, 
thus it has been fixed equal to O.Iujq- The effect of the damping can enter our calculation 
substituting in eq.(|58D Q^^^ for Q given by eq. (^5|) . 

Considering the two terms of eq.(|60D, the correlation function can be written as 



^aA^^n) = n« (^tu„) + ng),(zw,). 



(61) 
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The first term reads 



2.2 -2St I ( 'S'o + 1/2 



y)) (^)Esi^'(^«)^'e^""^(k,-r)^(k,r) (62) 



\V 2^ + 

and the second one is obtained retaining only the main autocorrelation term i = i' and 
6 = 6' 



'2f 



' So + 1/2 ' 
25 + 1 



1_\ J2 re--^(k,-r)^(ki,r)A(r) (63) 



where /i(r) is given by eq.(0). We stress that substituting Q^'^^ for Q is necessary to have a 
non vanishing Il'^\{iuJn)- We perform the analytic continuation icUn to + i6, and, clearly, 
the conductivity can be expressed as sum of two terms 



(64) 



The first term 'Sia^'^'^^ represents the band conductivity. Here the charge transfer is affected 
by the damping due to interactions with the lattice and spin fluctuations, but it is not ac- 
companied by processes changing the number of phonons. On the other hand, the second 
term ^a^"^"^^ in eq. (|64D derives from inelastic scattering processes of emission and ab- 
sorption of phonons and it is marked by the apex "incoherent" since the particles lose their 
phase coherence because of this phonon-assisted hopping. As in the spectral properties, we 
stress the appearance of two physically distinct contributions. 

In appendix we report the details of calculation for the optical conductivity. The band 
conductivity is derived as 



where C(^, u) is 



2S 



(66) 



and h{^) reads 



MO 



(67) 
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The latter term of the conductivity becomes 
(2eH 



UJ] 



//S'o + l/2V^ 



~2St 



-2St 



2S+1 



-oo r+oo 
J — oo 



h{z) 



'5o + l/2 



oo J — oo 
\ 2 



5o + l/2^ 



25+1 ;/ \\2S + l ) I 

/+00 r+oo 
/ d^ig{09{^i)[nF{^-uj)-nF{0]C{^,^uUj) 
-oo J —oo 

where g{^) is the density of states ([TgD and the function R{^,^i,u) is given by 

R{C, ei, ^) = 2 2 h{z) sinh j [JKe, 6, ^) + ^^Ke, 6, ^)] . 

We notice that Ji{^,^i,u!) 



(68) 



(69) 



Ji{^, 6, ^) = C'(e, ei, ^ + l^o) [ripi^ - too - ^) - npi^ - too)] [A^o(too) + upiO] (70) 
and Hi{^,^i,uj) 

HiiC, ^1, cj) = C(e, 6, ^ - too) [ni.(e + too - cu) - ni.(e + too)] [Noilcuo) + 1 - ^f(0] (71) 

describe phonon absorption and emission processes, respectively. 

In the limit of high temperatures (T > 0.39ci;o) and for small polaron excitations, the 
incoherent absorption is prevalent. In this case the conductivity consists of a sum of narrow 
Lorentzian functions centred on the points huq respectively.ii'0 We perform an expansion 
suitable for high temperatures and in appendix we derive the approximate expression of the 
conductivity 



( ff \ / ( So + i/r 

[ CO \ \ 2S + 1 , 




exp 



exp 



^ I I -\2 



p(l - p) (72) 



where p is the electron concentration. 

We have checked the internal consistency of our approach by means of the sum rule 



^ 



(73) 
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where Tq « is the component of the kinetic energy. The left-hand side of eq. (|73D is calculated 
by using the conductivity spectra, whereas the mean value < Tq, „ > can be evaluated 
by performing the two canonical transformations (@J^) and making the decoupling in the 
electron, phonon and spin terms through the introduction of Htest (§)• We get 

< T„,„ >= - / digmnpii) (74) 



3 

realizing that the two sides of eq.(|73D differ in a few per cent. 

In Fig. 6 we report the calculated conductivity spectra for x = 0.3 at different tempera- 
tures. In the large polaron ferromagnetic phase, a is given by eq.(Q), in the small polaron 
paramagnetic phase by eq.(|72|), in the regime of coexisting phases by the two preceding 
conductivities by means of the eq.(|l^). With rising temperature, a transfer of spectral 
weight from low to high energies takes place. The crossover energy can be estimated about 
Uc = 9uj(). When the temperature increases, the large polaron Drude term makes smaller 
and at T = Tc the optical response is characterized by the small polaron absorption band 
peaked approximatively around 2g'^ujQ. Thus changes in the quasi-particle excitations of 
the system are traceable in the optical response that, for T > T^, show the opening of an 
optical gap. At low temperatures the conductivity shows a clear Drude peak below ojq and 
a quiet asymmetric absorption band with a long tail above its peak position at higher ener- 
gies. These two types of optical response are due, respectively, to a coherent motion and an 
incoherent absorption band of large polarons. In the inset of Fig. 6 the Drude peak is shown 
. As T decreases, the coherent response increases continuously. We stress that the results 
shown in this figure are all consistent with experimental spectra.0Ei 

In Fig. 7 at T = O.lTc the optical conductivity is decomposed into two components: band 
(dotted line) and incoherent conductivity (dashed line). It is worth-while to notice that the 
band conductivity provides also the main contribution to the polaronic absorption band. In 
the inset of Fig. 7 the Drude term is fully shown. At low temperatures it can be obtained by 
means of the equation 

^^CyOrudeiuj) = A ^J[^] ^ (75) 
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where A is a proportionality constant and r(/i) is the scattering rate of the preceding section 
calculated at fi. This inset emphasizes the spectral weight of the absorption band too. 

In order to be more quantitative in the comparison with experimental data, we have 
evaluated the effective carrier number below the crossover energy Uc 



that is proportional to the spectral weight below this cutoff energy (see Fig. 8a). In eq.([76|) 
rrie is the free electron mass, e is the electron charge and Nmu is the number of Mn atoms per 
unit volume. As T increases, iVg// decreases assuming a nearly constant value for T > Tc. 
We stress that the calculated quantity (circles) agrees with the experimental effective carrier 
number deduced by ref.llllli and reported in the same figure (diamonds). We note that the 
agreement gets worse for T > Tc. Indeed, within this temperature range, in addition to the 
small polaron band, other terms, such as interband transitions between the exchange-split 
conduction bandsii, could contribute to the optical absorption. 
We have also estimated the quantity < u > (the first moment) 



as a function of the temperature (see Fig.Sb). It measures the characteristic energy of the 
polaron excitations.B Thus < u > assumes small values in the large polaron phase at low 
temperatures and increases gradually in the phase separation regime tending to the value 
g'^ujQ in the small polaron phase at high temperatures. 

Finally in Fig. 9 we plot the resistivity p = 1/a , where a results from the conductivity 
cr(co') in the limit uj 0. We obtain a further confirm of two different behaviors: metallic 
{dp/dT > 0) and insulating {dp/dT < 0). In logarithmic scale the behavior of activated 
resistivity^ for T > Tc is emphasized. 




(76) 



< UJ >= 



(77) 



duj^aa,ai^) 
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V. SUMMARY AND CONCLUSIONS 



We have discussed spectral and optical properties mainly for x = 0.3 and as functions of 
the temperature. The polaron formation and the small-to-large polaron crossover through 
the phase separation regime turn out to play a crucial role in understanding the properties 
of manganites. 

It has been stressed that the phase separation regime induces a robust pseudogap in 
the excitation spectrum of the system. Furthermore the pseudogap features show a direct 
relationship with the metal-insulator transition and with the crossover from the coherent 
large polaron dynamics to the incoherent small polaron one. Our results turn out compatible 
with the spectral weights found experimentally.0'0 

With regard to the infrared absorption, we have observed that, with decreasing T, our 
optical spectra are characterized by a transfer of spectral weight from high to low energies 
filling up the low frequency optical gap present in the high-temperature phase. 

In the metallic ferromagnetic phase, at low temperatures, the system shows two types 
of optical response: a narrow Drude term and a broad absorption band due respectively to 
the coherent and incoherent motion of large polarons. These results obtained at low T are 
consistent with experimental conductivity spectra.El'Elil The agreement improves when a 
temperature-independent part is extracted from the experimental optical spectra.iilii This 
part is considered to be due to the "background" interband transitions between the O 2p and 
Mn 3d band. The reduced optical conductivity obtained by subtracting the temperature- 
independent termii points out that the broad band at low T is very similar to that shown in 
this paper. We also stress that orbital degrees of freedom could be responsible of the effect 
of enhancing the incoherent absorption. E3Ei On the other hand, we observe that electronic 
correlations can change the spectral properties but are somehow obscured in the optical 
conductivity within the ferromagnetic metallic phase.0 

At high temperatures, the infrared absorption is due to the incoherent small polaron 
motion. Spectroscopic evidences of the small polaron formation have been reported in the 



23 



high-T paramagnetic phase. EJilj However for frequencies higher than the peak of small po- 
laron absorption band, the effects due to the exchange-split bands and to the local Coulomb 
repulsion cannot be neglected.^ Indeed the experimental data can be fitted reasonably well 
with the small polaron band below 0.8 6^.0 Thus we can consider our results meaningful 
for frequencies up to the absorption peak. 
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APPENDIX 

In this appendix we report the calculations leading to the expressions of the optical 
conductivity 

^a^A^) = 3fJa2r'n^) + 3^<r'H^)- (78) 

This decomposition originates from the current-current correlation function that in Matsub- 
ara frequencies can be written as 

n„,„(zcu„) = n^i^j^^n) + ng)jia;„) (79) 

where the first term is 

n« (^-n) = 4eVe--- ^ (^1^) )' (^) E ^-'(^") r -)^(k, (80) 

and the second one reads 

with /i(t) given by eq. (|3^) . Adopting the Lehmann representation, the spectral function 
A(k, u) can be introduced, therefore, making the analytic continuation iun ^ uj + i6, the 
band conductivity becomes 
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-jX:sm^(M^(k) 



(82) 



where D{k) = L)(k, k) with 

D(k,kO = 



+00 g^l^' _ 

- — A(k, ci;')y4(ki, a;' — ci;)[nir(a;' — 00) — npioo')]- 
-00 47r 



(83) 



We also derive the latter term of the conductivity 



^n(y(a;) 



UJ 



2S+1 



UJ 



h{z) 



25+1 



25+1 



]V^)E^(k,ki) 



k,ki 



(84) 



where 



+00 



' l3uJol 



f{k, ki) = 2 5] Ii{z) sinh ^ [J^k, ki, a;) + Hi{k, ki, u;)] 



(85) 



with 



too (j^ijj' _ 

J;(k, ki,a;)= / - — 74(k, a;')74(ki, a;' — Zcx)o — a;)[ni7(a;' — Z(x)o — a;) — ni^(a;' — ZcDq)] >< 
-00 47r 



(86) 



and 



r+00 ^ 

if;(k, ki, a;) = / - — A(k, cu')A(ki, a;' + ZtDo — '^)[^f('^' + ^'^o — <^) — ?T'f((^' + ^'^o)] >< 

-00 47r 



X [Aro(too) + 1 - nj.(u;')] 



(87) 



We observe that J/(k, ki,a;) and i7/(k, ki,a;) describe phonon absorption and emission pro- 
cesses, respectively. Using the convolution product of Lorentzian functions, the band con- 
ductivity becomes 
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^^ta'\-^)=[^)e-''^{[^^^ (88) 
where C{^,uj) = C{^,C,,uj) with 



and 



M0 = (^)Esi^'(M5(^-a). (90) 



Since ^'{y), the Fourier transform of h{^), is calculated 



we can express h{^) in the following manner 

2nteffJ~^ \2teffJ y 

On the other hand, the latter part of the conductivity reads 



2eH 



2+2' 



-2St 



'Sq + 1/2V\ //5o + 1/2\\'' 



u J \ ■ \\2S+l ) I \\2S + l ) I 

/+00 r+oD 
/ d^ig{09i^i)[nF{^ - ^) - n^(0]C(e,ei,c^) (93) 
-oo J — oo 

where g{^) is the density of states (|16D and the function R{^,^i,u) is given by 

R{t ei, c^) = 2 E ^K^) sinh (^^^ [ J,(e, ei, ^) + ei, ^)] (94) 

with 



M^, 6, = C'(e, ^1, ^ + /^o) [nri^ - too - a;) - - too)] [A^o(too) + (95) 



and 



Hi{^, Ci, = 6, ^ - too) + too -co)- ripi^ + luo)] [iVo(too) + 1 - ripiO] ■ (96) 
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In the limit of high temperatures, we concentrate on /i('r) of eq.(|8TD carrying out the 
following approximation 



e ^'^'^ exp \ z cosh 



exp(— 2;rj exp I — r 



where 



(97) 



(98) 



This expansion is valid as long as (^yf ) ('^) < 1; hence T > O.SQcuo. The response 
of the small polaron excitations is characterized by the incoherent behavior, therefore the 
conductivity becomes 



eH^\ I (S^ + l/2\ 



iTlp f 1 



u ) \ \2S + l J I 
X [1 — nF{uj")] nir(a;')< exp 



E 

k,ki 



27 l-oo ^^(k,^ )^(ki,^") X 



^{UJ + uj' - uj" - zf 
Az 



exp 



/9 



Az 



(99) 



We consider ~ —g^ujQ — ;U obtaining 



(eH^\ //5o + l/2' 



hl3 

— lexp 

z 



\\ 25+1 / 
where p is the electron concentration. 



exp 



p(l-p) (100) 



FIGURE CAPTIONS 

Fl (a) The phase diagram corresponding to t = 2ujq, g = 2.8 and e = O.Olt. PI means 
Paramagnetic Insulator, FM Ferromagnetic Metal and AFI AntiFerromagnetic Insu- 
lator. The areas PI + FM and AFI + FM indicate regions where localized {PI and 
AFI) and delocalized (FM) phases coexist, (b) The fractions of volume iVi/V)ferro 
and (y2/V)para at X = 0.3 as a function of the temperature (T is expressed in units of 
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F2 The spectral function (corresponding to t = 2u!o and g = 2.8) at x = 0.3 as a function 
of the energy (in units of t^o) for two different temperatures: (sohd hne) T = 0.05 and 
= 0, (dashed hne) T = 0.47 (T is expressed in units of luq). 

F3 The renormahzed density of states at x = 0.3 as a function of the energy (in units of 
Uq) for different temperatures. In the inset the renormahzed density of states at /i as 
a function of the temperature. The resuhs are obtained for t = 2ujq and g = 2.8. 

F4 The renormahzed density of states at T = 0.390^0 as a function of the energy (in units 
of a;o) for different densities. In the inset the renormahzed density of states in a small 
range of energy around the chemical potential /x. 

F5 (a) The scattering rate (corresponding to t = 2ujo and g = 2.8) at x = 0.3 as a function 
of the energy (in units of u^o) for different temperatures, (b) The scattering rate (solid 
line), 1-phonon rate (dotted line), spin- fluctuation rate (dashed line) and multi-phonon 
rate (dotted-dashed line) at x = 0.3 and T = 0.65Tc as a function of the energy (in 
units of Uq). 

F6 The conductivity up to 18 c^o at different temperatures. In the inset the conductivity 

2 

up to 1 uq at different temperatures. The conductivities are expressed in units of 
with c hole concentration and = ^, and they are obtained for t = 2uJo and g = 2.8. 

F7 The conductivity (solid line), the band conductivity (dotted line) and the "incoherent" 

2 

conductivity (dashed line) at T = O.lTc in units of where c is the hole concentra- 
tion and ^ = In the inset the conductivity (solid line) and the Drude term (dotted 
line) at T = O.IT^ in the same units of the figure. 

F8 (a) The effective carrier number (circles) below the crossover energy uJc compared to 
the experimental effective carrier number deduced by ref.0Jii as a function of the 
temperature, (b) The first moment (in units of Uq) as a function of the temperature. 
The results are obtained for t = 2u!o and g = 2.8, in the case (a) specifying ujq = 50 
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meV and the lattice constant a — 0.4 nm. 

F9 The resistivity as a function of the temperature (we have used cuq — 50 meV and the 
lattice constant a — 0.4 nm). 
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